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ABSTRACT 

A  general  class  of  estimators  is  developed  for  improving 
upon  best  scale  invariant  estimators  of  two  or  more  arbitrary 
scale  parameters  (or  powers  thereof)  for  arbitrary  positive 
distributions  with  sufficient  moments  under  weighted  squared 
error  loss  function.  The  technique  i6  to  compute  the  risk 
difference  in  terms  of  moments  of  the  distribution.  Some  condi¬ 
tions  are  obtained  under  which  the  maximum  improvement  i6 
possible,  and  the  form  of  the  estimator  can  be  chosen  to  achieve 
this  maximum  along  any  specified  ray. 

The  result  is  then  extended  to  the  estimation  of  a  linear 
transform  of  the  parameter  vector.  Finally,  some  examples  are 
given  with  numerical  calculations  to  obtain  the  amount  of  risk 
improvement. 


1.  INTRODUCTION 

There  has  been  considerable  interest  in  improving  upon  the 
standard  estimators  in  multivariate  estimation  problems.  Since 
the  celebrated  work  of  Stein  (1955),  numerous  results  have  been 
proved  which  show  that  the  presence  of  Stein  effect  is  just  part 


■M»wro«i«PHnia«»po»»«inBr«H»nEnHnm»iargii  E»nnTX>iuxTXX*vinn>grsinaKa»gain»uwuwgwtt»nwH»uwCTCTrowgir»< 
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of  a  very  general  phenomenon,  which  has  little  to  do  with  the 
exact  form  of  the  loss  function  or  the  underlying  distribution. 

The  simultaneous  estimation  of  location  parameters  was 
investigated  by  several  authors  during  the  past  decade  culminating 
in  a  unified  discussion  by  Shinozaki  (1984).  The  improved  estima¬ 
tion  of  scale  parameters  from  exponential  families  has  also  been 
studied  recently.  The  major  results  in  this  direction  are  obtain¬ 
ed  by  Hudson  (1978),  Berger  (1980),  and  Ghosh,  Hwang  and  Tsui 
(1984).  The  basic  technique  for  obtaining  improved  estimators  of 
scale  parameters  in  exponential  families  is  to  obtain  an  unbiased 
estimator  of  the  risk,  using  integration  by  parts  technique  and 
solving  a  differential  inequality.  Recently  DasGupta  (1986), 
however,  obtained  improved  estimators  of  gamma  scale  parameters 
without  a  variational  argument  by  proposing  an  estimator  which  i 6 
a  function  of  the  geometric  mean . 

In  this  paper  we  extend  DasGupta  (1984,  1986)  in  several 
directions.  We  consider  the  problem  of  simultaneous  estimation  of 
arbitrary  powers  of  scale  parameters  from  arbitrary  independent 
positive  valued  distributions  under  weighted  squared  error  loss. 
All  that  i6  required  is  that  moments  of  a  specified  order  be 
finite . 

Suppose  X  *  (X^,...,X  ),  p  ^  2,  where  X^  has  density  fg  (xj 

1  P  i 
m  t—  f .(x./6 .) ,  x.  >  0,  6.  >  0,  i  ■  l,...,p.  Also  suppose 
6.1111  1 


*  2s . 

Efl  (X.  l)  < 
.  1 
1 


e 


(e 


8i 


*  •  I 


i  *  l,...,p.  Our  problem  is  to  estimate 

s 

6  **)  under  the  loss  given  as 


L(a,6s) 


pc  s.2 

I  w.6.  (a.  -  6.  )  . 
i-1  1  1  1  x 


(1.1) 


It  can  be  easily  checked  that,  under  the  loss  (1.1),  the  best 

invariant  estimator  of  68  is  6®(X)  *  (6^(X) . 6®(S)),  idiere  6?(X) 

1  P  i 


is  given  as 


•  •  d.  o  • 

where  r.  „  ■  E(X.X)/E(X.  ),  the  expectation  E(*)  taken  under 

1.8. .  1  1 

i  i 

6,  ■  1.  In  Section  2,  we  propose  a  class  of  improved  estimators 
1  0 

which  dominate  6  (X).  Another  class  of  estimators  are  obtained 
which  incorporates  a  prior  guess  about  the  parameter  vector. 

In  Section  3,  we  study  the  component  risk  behavior  of  certain 
"best"  estimators  in  the  class  obtained  in  Section  1.  By 
aggregating  the  component  risk,  we  can  obtain  an  explicit  form  for 
the  improvement  in  risk  and  identify  the  direction  in  which 
maximum  improvement  is  attained. 

Section  4  is  devoted  to  the  estimation  of  a  linear  combina¬ 
tion  Jt'e8  of  the  parameter  vector,  where  the  vector  i  . . .  t  )' 

g  P 

is  p  x  1,  Several  estimators  are  obtained  which  dominate  aX 
where  a  *  (a.,,,. a  )'  is  p  x  1  vector  and  X8  ■  (X  , . . .  ,X  **).  Some 

1  p  i  p 

examples  are  given  which  include  estimation  of  £6.. 

i  1 

Finally,  in  Section  5  we  give  examples  of  simultaneous  esti¬ 
mation  of  ratio  of  independent  normal  variances  and  estimation  of 
ranges  of  rectangular  distributions.  Some  numerical  studies  are 
performed  which  give  the  percentage  improvements  in  risks. 

2.  IMPROVED  ESTIMATORS  OF  6* 

g 

In  this  section  we  will  develop  estimators  of  6  which  will 
dominate  the  best  invariant  estimator  (1.2)  under  the  loss  given  in 
(1.1)  for  p  ^  2. 

Brown  (1966)  and  Brown  and  Fox  (1974)  provide  perhaps  the  most 
general  conditions  for  admissibility  of  the  best  invariant  scale 
parameter  in  the  one  dimensional  case.  By  transforming  the  scale 
problem  to  the  equivalent  location  problem,  it  is  straightforward 
to  verify  that  r  .X  is  an  admissible  estimation  for  6  under 

6  }  Z8 

_  2 

squared  error  loss  if  EX  vtnX) < «° where  the  expectation  is  taken  at 
0*1.  That  is,  referring  to  conditions  (a)  -  (d)  on  p.  808  of 
Brown  and  Fox,  since  the  best  invariant  estimator  is  unique  (a)  and 
(b)  are  apparent.  The  expectation  condition  immediately  implies 


(c)  and  (d)  as  well  by  appealing  to  Lemma  2.3.3,  p.  1105,  of 
Brown  (1966). 

Now  we  will  develop  the  following  notation: 

(1)  m.  «  E  X“; 

l  ,a  lx 

(2)  ri,a..e.  "  Bi,a./,ni,6.; 

(3)  ga(X)  -  llX^-1;  (2>1) 

(4)  E.g  (X)  -  Urn.  -  p  ; 

la  *1.0.  a 

J  J  J 

(5)  Elt'<»  -  n„  -  P 

J  J 

With  the  above  notation,  it  follows  that  E.  (xT)  *  eTm.  , 

6i  i  i  i,r 

i  *  l,...,p.  Classes  of  improved  estimators  are  given  in  Theorem 
2.1,  but  first  we  note 

Lemma  2.1.  For  any  a ^ ,  j  »  l,...,p, 

a./p  a 

Pne.J  /E6.J  <  1. 

j  J  j  J  " 

Proof.  The  arithmetic  mean-geometric  mean  inequality. 

Lemma  2.2.  Under  one  of  the  following  conditions, 

(i)  c  >  0  and  t  >.  0, 

(ii)  c  <  0  and  t  <  0, 

(iii)  c  *  0  and  arbitrary  t, 
the  following  inequality  holds: 

ne!/p  re’r/ze':*6  <  1. 

3  J  j  J  j  J  " 

Proof.  The  proof  for  (i)  and  (ii)  follows  from  Holder's 
inequality. 

Lege  a  2.3.  r.  ..  as  defined  in  (2.1)  is  an  increasing  func- 
tion  of  s  if  t  >  0,  a  decreasing  function  of  s  if  t  <  0. 

Proof.  It  is  clear  that  X.  is  monotone  in  X.  and 
■  l  i 

h.(x,t)  *  x.f .(x)//x.f .(x)  dx  is  MLR  in  X.. 

1  1  x  '  x  x  i 
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is 


t+8  t 

Thus  EX^  /EX  ^  is  increasing  in  s  if  t  >  0,  and  decreasing 


in  s  if  t  <  0. 

We  are  now  ready  for 

Theorem  2.1.  Consider  the  estimator  6(X)  given  componentwise 


as 


6.(X)  -  6?(X)  -  bsgn[(ai  ♦  B.  -  si)si)Xilga(X) ,  (2.2) 


i  *  1 ,  —  ,p(>^  2)  , 


where  cu  ♦  B^  -  s^  j4  0  for  all  i  *  l,...,p  Then  6(X)  dominates 


d^tX)  in  terms  of  risk  if  either  (1)  or  (!')  and  (2)  hold: 


(1)  a.  *  (c.  +  B.  +  s.)/p  and  B.  *  -c./2,  i  ■  l,...,p;  (2.3a) 

l  l  l  i  i  i  r 


(1')  All  c^  equal,  all  B^  equal,  all  s^  equal  and 


*  a  *  (s-B)/p; 


and 


(2)  0  <  b  <  2pad(1)/p2od<2) 


(2.3b) 

(2.4) 


where  d^^  and  d^2^  are  positive  and  defined  below. 


Proof.  Let  4(6)  ■  R(6,d)  -  R(8,dg),  the  risk  difference.  It 


is  sufficient  to  show  that  4(6)  <  0.  Now,  using  the  loss  (1.1),  it 
follows  that 


c.  ,  2B.  , 

i,,2„  i  2 


s,  6 


4(6)  «  E.Uw.e.^btX.  A  g‘(X)  ♦  2r .  ,  X.1  b.X.1  g„ 

0,11  ii  a  i,8. ,28.1  ii  a 


(x) 


s.  B . 


2b, 


ieilxilga(X)}],  b.  -  -b  sgn[(a.  ♦  6.  -  sJs.], 


c.  •  28.  p 

*  Zw.6.  {b.6.  g  (6)(rim.  -  )m.  «r  _  v/m.  - 

•  ii  ii  ®2a  .  ,j,2a.  i,2(a.+B.)  i,2a. 

l  l  J  j  i  l  i 


B.+s . 

+  2r.  _  6 . 1  S.g  (6)(IIm.  )m.  .  /m. 

i,s. ,28.  i  l  a  .  j ,a .  l, a. ♦B.+s.  i,a. 
’i  i  J  J  ill  i 


B.  s. 

2b. 6. 1  g  (6)(Ilm.  )m.  _  /m.  } 

ii  a  .  j,a.  i,a.+B.  i,a. 

j  j  ii  l 


6 


iBWiWCTBBBBPaPBCiPPnqwnqngnqrawqnannEiatCTnfraTiErainnnniPTgwCT^ 


^Wiei1{bi6i  1*2a(e)p2ari,2(a.+e.),2a. 

1  ill 

6. +8. 

♦  2r .  ,  b.e.1  g  (®)P  r.  0 

i,6. ,2s.  ii  °a  a  i.a .+B.+8 . ,a . 
i  i  1111 

6. +e . 

-  2b. 6.  xg  (0)p  r  ,0  } 

i  i  ®o  a  a.+B. .a. 

1  X  I 

2 

*2c,(6)c2!.  1  "iVi  ri,2(a.+B . )  ,2a . 

1  111 

C.  +0.+S. 

+  2g  (e)p  z  w.e.1  1  1  b.d. 

a  a  ^  ii  ii 


(2.5) 


where  d . 


ri,s.,2s.  ri,a  .+6 .+s . ,a .  ri,a.+B.,a.’  *  1*---*P* 

i  i  1111  ill 


Now,  using  Lemma  2.3,  it  follows  that  for  b.  >  0  if  a.  + 

<  s.,  then  d.  <0;  and  if  a.  +  B.  >  6.,  then  d.  >0.  Similarly 

l  i  ill  l  ' 

for  s.  <0,  if  a.  +B.  <0,  then  d.  <0;  and  if  a.  +  B.  >  ».,  then 

ill  l  ill 

d^  <  0.  Note  that  b.d^  is  always  negative.  Now,  define 

d^  *  min  w.  Id.  I 
i  11 

4  ri,2(a,.B.),2».- 

1  1  1  1 


It  follows  from  (2.5)  that 

.2.(2), 


c.+2B. 


c.+B.+s. 

d>™  ill 


A(6)  <  g,„(e)p  b  dw;ie 1  ~  2g  (6)p  bdvuie. 

2q  2a  i  °a  a  l 


-  g  (e)ie.1  1  1{b2p,  d(2)K(e)  -  2bd(1)P  ), 

a  i  Za  a 


(2.6) 


id»ere 


c.+2B.  c.+B.+s. 

K(e)  -  g  (e)re.1  Vie.1  1  V 


If  (1)  of  (2.3a)holdB,  then  from  Lemma  2.1  we  have  K(6)  <  1, 
i.e..  Me)  <  0  if  (2.4)  holds.  If  (1’)  of  (2.3b)  holds,  then  with 
Lemma  2.2  we  take  t  ■  s  -  B  and  a  *  (s-6)/p.  Note  that  a  -  B  ¥  s, 
i.e.,  ( 8—b) / P  ♦  B  f4  8  requires  B  8  snd 


if  c  ♦  28  >  0  ^  0,  i.e. ,  6<  a  -c  <  28  <  2s,  if  c  >  -  2s 

if  c  + 26  <  0  — >  i  <  0,  i.e.,  6  >  •  “>  2a  <  26  <  -c,  if  c  <  -  2s 
c  ♦  28  *  0  c  -  2s. 

Again  from  (2.6),  it  followa  K(0)  <  1  ao  that  4(0)  <  0  if  (2.4) 
holda.  Thia  completea  the  proof  of  the  theorem. 

Remark  2.1.  Theorem  2.1  ia  a  more  general  veraion  of  Theorem 
2.1  in  Das  Gupta  (1984). 

Remark  2.2.  From  the  proof  of  Theorem  2.1,  it  follows  that  if 

for  any  i,  a.  ♦  8.  ■  s.,  i.e.,  c.  ■  -2s  ,  the  dominance  result  does 
J  1  1  1  x  £ 

not  follow.  Thus,  under  the  invariant  quadratic  loss,  our  result 

does  not  quite  give  an  improved  estimator.  However,  if  a  minor 

®  i 

restriction  is  taken  on  the  parameter  space,  say  min  6.  >  c,  where 

i  1 

e  is  a  preassigned  positive  number,  then  estimators  of  our  form 
will  dominate  the  best  invariant  estimator  under  the  invariant 
quadratic  loss. 

Now  we  will  consider  further  generalization  of  Theorem  2.1. 

For  example,  suppose  we  have  a  prior  guess  for  6,  say, 
eQ  *  (0Q1 , .. . «eop) *  How  “ighf  we  incorporate  this  information  in 
our  improved  estimator?  The  following  theorem  gives  the  improved 
estimator  in  this  case. 

0Q 

Theorem  2.2.  Consider  the  estimator  6  (X)  given  component¬ 

wise  as 

0  0 
6.°(X)  «  6?(X)  ~  b  sgn[(a  +B.-6.)s.] (X./0-.)  1  g  (X/6-), 

x  i  i  x  x  i  i  0i  0 


i  m  1 » •  >  •  ,p. 


(2.7) 


Then  6  ®(X)  dominates  6®(X)  in  terms  of  risk  if  either  (2.3a)  or 
(2.3b)  holds  and 


(2.3b;  holds  and 

« < » <  vvv  S.*  • 

where  d  and  d  are  positive  and  given  below. 


(2.8) 


pi 

Proof.  Define  k^  ■  6q£  ,  i  *  l,...,p.  Now,  we  can  write 


;d 


si 


$} 


.*tl 


6 .  (X)  -  r.  _  X.  +  b.k.g  (6n)g  (X)X.‘ 
x  1 ,  s . ,  2  s .  i  i  i°o  0  “a  l 


where 


~  b  sgnKaj+B^s^ij]  ,  i  - 


Nov,  it  clearly  follows  that  the  risk  difference 


^  •  9  9  9  9 


Me)  -  z  Vi  * biE-a<eo )kiExi  ‘(g^x)) 


gAOjk.Eiv'^g  (X)) 


+  2b. r.  ,  g_C6.)k.E(v.  X.  g 

l  i,s^,2sj  -a  0  l  l  i  a 


"  2bieIV-a(e0)E(XiVX)>} 


c.+2B.  , 

,  i  i  v2, 


*2a<e)*-2a(80>C2ai”i6i  bikiri,2(a .*6. ) .2« . 

1  111 


+  28a(6)pa8-a(e0^  Vi 


c.+B.+s. 

ill 


bi^ri,s.,2s.  ri,o.+6.+8. ,a .  ri,a .+B . ,a . ^ * 
i  l  l  i  i  i  l  i  i 


(2.9) 


Define 


d  *  min  v.k.  d. 

l  l 1  l 1 
i 


*(2) 

d  ■  max  v.k.r.  v  _ 

l  i  i,2(a.+0.),2a. 
l  ill 


Now  following  calculations  similar  to  that  in  proof  of  Theorem  2.1, 


it  follows  that 


9  h( 9  ^  C .+2  B . 

Me)  <  B2a(e/e0)p2ab2d  {  }  zei1 


,r.x  c.+B.+e. 

-  28a(e/60)pad  E6il  1  * 

l 


(2.10 


As  in  the  proof  of  Theorem  2.1,  it  follows  that  if  (2.8)  holds, 
then  A(6)  <  0.  This  completes  the  proof. 


V.’ 


Remark  2.3.  Das  Gupta  (1984)  considers  a  special  case  of 
estimators  which  in  our  notation  would  extend  (2.2)  to 


» . 

6 . (X)  -  6T(X)  +  b  agn[(a.+B.-e.)s.)X.lg  (X)6  * 


i  l  l  l  l  a 


(2.11) 


It  is  clear  that  by  transformation  (2.11)  is  equivalent  to  (2.7) 


Remark  2.4.  Notice  that  if  6^  ■  1,  i  ■  l,...,p,  then  the 


estimator  (2.7)  is  the  same  as  that  of  (2.2). 


Remark  2.5.  It  is  noteworthy  that  the  range  of  b  in  (2.5) 


depends  also  on  the  values  of  6^. 


Remark  2.6.  The  dominance  result  holds  even  if  X.'s  arise 
—  —  i 


from  different  families  of  distributions,  say,  one  from  exponential 


and  the  other  from  uniform. 


3.  COMPONENT  RISK  AND  MAXIMUM  IMPROVEMENTS 


In  this  section  we  study  the  component  risk  behavior  of  the 


estimators  given  in  (2.2)  and  in  (2.7)  when  f^(x)  *  f(x), 


i  *  1 . .  In  this  case,  we  take  r.  e  ■  r  and  m. 

r  i,a,B  o,B  i,8 


The  following  theorem  gives  the  maximum  improvement  of  the 
estimator  of  the  form  (2.2)  when  b  is  the  midpoint  of  (2.4)  and  the 
direction  at  which  the  maximum  attains. 


Theorem  3.1.  Suppose  X.  is  distributed  with  pdf  fg  (x^) 

i 


■  •g—  f(x/6.),  i  *  1 . p.  Consider  the  estimator  defined  component- 

.  i 

wise  as 


6 .  (X)  -  r  ,  Xs  -  bX?  H  X“,  i  «=  l,...,p, 
i  s,2s  1  1  .  j 


(3.1) 


With  b  “  ra,2ad/r2(a+8),2a 


where 


d  ■  r 


s,2s  ra+8+s,a  ra+8,a‘ 


(3.2) 


Then  the  maximum  improvement  in  risk  when  the  loss  is  (1.1)  with 


c.  ■  c  and  w.  ■  1  is  attained  when  6.’s  are  equal,  say,  to  6  and  the 
l  l  i 


•mount  of  improvement  is 


2  2p  .c+28. 

r2«,2{.*e)“« PB  /bV 


(3.3) 


Proof .  The  difference  in  risk  for  the  ith  component  is  A^(6) 


where 


A.(e)  -  R(6?,ef )  -  R(6.,e®) 

i  ii  ii 


-e^tE{b2x?6nx^a  -  2bx?  nx°(r  ,  xe.  -  e®> >] 

i  i  .  j  i  .  j  s,2s  i  i 

j  j  J 

cf.2  2(a+B),  (ah2  p  , 

6.[b  6.  (n-i  >  V2(a+B)/o2a 


-2b0?+a+8(r/<?^  )mP(m  .  r  ,  /m  -  m  /m  )] 
l  -l  a  a+B+s  s,2s  a  a+B  a 

(a) 


n  e“) 


(where  n 

-i  ...  J 

J**- 

-RCrh2fl2(a+e)<’n(0,^2mP  r 

-e.[b  e.  (n_.  )  m2a  *2(a+e)>2a 


-2bea+6+8(r,<(?))  mPd] 


-1 


(3.4) 


Substituting  b  in  (3.4),  we  obtain 


A.(6) 


„cf,  p  ...  .o+B+b  ,  p  (a) 

6 .  2rr  .  (d/r_,  ,0x  ,  )  6.  dmrn  . 
i  a, 2a  2(a+B),2a  i  a  -i 


/  2p  ,2 

Va,2a/r2(a+B) 


2<n<“))' 

“I 


2(a+6)  p 

i  m2ar2(a+B),2a 


1 


e?+2s  d2  —  {e6+a-en(a)(2  -  e?+ot"s  n(a))>.  (3.5) 

1  r2(a+B) ,2a  1-1  i  -i 


From  (3.5)  it  follows  that  the  maximum  improvement  occurs  when 

ea+B-s  n(«)  .  ^ 
l  -i 

Since  this  is  true  for  each  i  *  l,...,p,  maximum  overall  improvement 
occurs  when  all  6^  are  equal,  say  to  6,  and  the  total  improvement  by 
summing  (3.5)  over  i  is 


11 


4S.,2(»»)  "a1"  »e 


2P  (3.6) 


This  completes  the  proof  of  the  theorem. 

We  note  that  stunning  (3.5)  over  i  yields  the  total  improve¬ 
ment  of  6  for  general  6.  The  component  risk  associated  with  6  is 
clearly  0?  8(1  -  m^/n^)  vh«ice  the  total  risk  for  is 

m2 

(i  -  !i_  )  re‘+28  .  (3.7) 

"2s 

If  we  define  the  percent  relative  improvement  in  risk  as 

R(6,6  )  -  R(0 ,6) 

PI  -  . - = -  x  100  ,  (3.8) 

R(0,6O) 

then  in  the  case  where  all  6^  are  equal  (3.8)  does  not  depend  on 
6  and  becomes 

t2p 

(3.9) 


d2  ®2a 


(m  Y 
a 


1  - 


ms  "2(0+8)  ("2o) 


m. 


2  s 


We  recall  from  (2.3b)  that  in  (3.1)  a  *  (s-0)/p.  Consider 
the  limiting  PI  86  p  -*■  ®.  If  EXC^  <  «  for  some  t  >  0,  this 


insures  that  lim  (m„  )  P  m2p 

2a  a 

p  00 

becomes 


Hence  the  limit  of  (3.9) 


(m„  r  .  -  m  ) 

B+S  8,28 

1  -  m  r  „ 
s  8,2s 


(3.10) 


In  the  case  when  B  ■  0,  s  ■  1,  this  reduces  to  1  -  m^  r^  2 
as  in  Proposition  4.1  of  Das  Gupta  (1984). 

Remark  3.1.  Theorem  3.1  is  a  more  general  version  of 
Das  Gupta  (1984)  Proposition  3.1. 


Theorem  3.1  obtains  maximum  improvement  along  the  ray 
(6 ,6,...  ,6)  =  61.  Suppose  we  have  prior  information  leading  us 
to  aeek  maximum  improvement  along  the  ray  defined  by  a  specified 
vector  6q.  If  we  use  the  estimators  in  (2.7)  or  (2.11),  we  will 

find  that  maximum  improvement  occurs  again  along  the  ray  61. 

Rather,  we  extend  Theorem  3.1  using  estimators  slightly  more 
general  than  those  in  (2.7).  Ue  retain  the  same  setup  as  in 
Theorem  3.1  and  consider  estimators  given  componentwise  as 

6i(X)  *  rs,2s  Xi  "  bi(Xi/0Oi>B  n(Wa*  1  *  1 . .  (3,11) 

In  order  that  such  estimators  dominate  6^,  we  need  the  risk 
difference 

«-2»<V  I2(a.e).2<,»L£bier26 

*  V9)  «-a(V  "a  d  1  vr6“  ‘  0 

with  d  as  in  (3.2).  Let  b.  *  h(6A.)b  where  h(*)  >  0  and  let 

h  *  min  h(0^),  h  *  max  h(6^  J.  Mimicking  the  argument  of  Theorem 

2.1,  if  b  is  between  0  and  2h  r£  r2(a+B)  2a’  then  ^.11) 

will  dominate  6^. 

Turning  to  component  risk,  if  we  take 

b  *  h  r**  „  d/(h)^  r.,  .  ,  then  the  difference  in  risk  for  the 

-  a,Za  Zva*B/ ,Za 

ith  component  is,  analogous  to  (3.5), 


30.2^2  r»,.2.  ,h(6  >i- 

1  r2(a+6) ,2a  °‘  <h)2  *  '*  °*  (h)2  1 


(3.12) 

The  maximum  improvement  occurs  when  h(0Q£)8®  *  is  constant  for 
all  i.  If  h(0Q.)  *  ^oi^  ’  occur*  slooR  the  ray  defined  by  ©q- 
We  summarise  the  above  as 


k*  WJFUWWVWWUlIWUWlWUqtlMIlWHJl \lfX  fuirjiiiiRH’iiK  CmrV  WWDTOWW*?  «7/«',vi'V’.Tfin  B 


Theorem  3.2.  Suppose  X.  distributed  with  p.d.f.  f_.(X.) 

_1 -  1  0X1 

.  £(X^/6.).  Let  0rt  be  an  arbitrary  positive  vector. 


x  l 


Estimators  of  the  form  (3. 11)  uniformly  improve  upon  6^  for  b. 


satisfying  conditions  given  below  (3.11).  Componentwise  risk 
improvement  is  given  by  (3.12)  and  is  maximised  along  the  ray  defined 


by  e0* 


4.  ESTIMATION  OF  A  LINEAR  COMBINATION 
In  this  section  we  will  consider  estimation  of  a  linear  com¬ 


bination  of  powers  of  scale  parametric.  Suppose  the  parameter 


function  to  be  estimated  is  y(.Q)  *  l  I.0.1  (we  take  £.  i  0  w.l.o.g.) 

.11  x 

X 


L(y(6) ,a)  *  (a  -  y(0))  , 


and 


6„(X> 


s . 

I  a.X.J 

j 


(4.1) 

(4.2) 


is  an  arbitrary  linear  estimator  of  y(6). 

The  following  theorem  gives  an  estimator  of  y(6)  which  ha6 


smaller  risk  than  6q(X)  under  the  loss  (4.1). 


Theorem  4.1.  Consider  the  estimator 


s./p 

6(X)  -  6n(X)  +  b  n  X.J  . 

0  i  J 


(4.3) 


Suppose 


and 


dj  “  £j  *jrj.8j(l+  l/p),s./p 


»  j  “  1 » •  •  •  »P  • 


'(1) 


min  d.  ,  d,  .  *  max  d.  . 

j  3  <p>  3  J 


Then  6(X)  dominates  6q(X)  under  the  loss  (4.1)  for  p  >_  2  if  one  of 
the  following  conditions  hold: 


(1)  dj  >  0,  Vj  *  l,...,p,  and  0  <  b  <  2Pg/p  Pd(1)/p28/p 

(2)  dj  <  0,  Vj  -  l,...,p,  and  2 pg/p  Pd(p)/P2e/p  <  b  <  0 
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If  any  d^  *  0,  j  *  this  approach  fails  to  provide  an 

improved  estimator. 

Proof.  Let  A(y(6))  -  R(y(6),6(X))  -  R(y(6) ,6Q(X) )  be  the 
risk  difference.  Then  it  clearly  follows  that 

9  2s. /p  a./p 

a(y(6))  -  bV  ,(iie.  J  )  -  2bY(e)p  ,  (n  e.J  ) 

Zs/p  j  j  s/p  j  J 

a.  s./p 

+  2bp  ,  l  a. 6.  (11  6.  )r.  .  , 

•/P  j  J  J  j  J  J.®j+6j/p.Sj/P 

2  2s. /p  s./p  s 

*  bp,  .  ne.J  -  2bp  ,  n  e.J  i  e.Jd. 

2s/p  j  J  ®/P  j  J  j  J  J 

Now  the  proof  follows  similarly  to  Theorem  2.1  using  Lemma  (2.1). 

Remark  4.1.  If  we  write  a.  ■  £.r.  _  ,  we  are  improving  upon 

J  J  J»u»sj 

the  corresponding  linear  unbiased  estimator.  Here 
d.  *  £.(1  -  r.  „  _  r.  _  »  _  ,_). 


l .  *  n.U  -  r .  _  r.  .  ,  ,  .  ) 

J  J  J,0,6.  J jS.+S./p,S./p 

J  J  J  J 


But  by  Lemma  2.3,  the  term  in  parenthesis  is  always  <  0;  »*»ence  if 

all  £j  have  the  6ame  sign,  the  unbiased  estimator  is  inadmissible. 

If  we  write  a.  *  £.r.  -  ,  we  are  improving  upon  the  corres- 

J  J 

ponding  linear  combination  of  best  invariant  estimators.  Here 

d.  *  £.(1  -  r.  _  r.  .  ,  .).  Now  by  Lemma  2.3  the  term 

J  J  J,Sj,2Sj  J,8>Sj/p,S  j/p 

in  parenthesis  i6  always  >  0.  Whence,  if  all  £j  have  the  same 
sign,  this  estimator  is  inadmissible  as  well. 

Remark  4.2.  Interestingly,  if  for  example  s^  *  1  and  all  £^  are 
equal,  say,  £  >  0,  i.e.,  we  are  estimating  £10,  then  with 

as  in  (4.2),  we  need  an  expander  if  £  is  sufficiently  large; 
whereas,  we  need  a  shrinker  if  £  is  sufficiently  small.  This 
paradoxical  dependence  upon  an  arbitrary  £  is  intuitively  reason¬ 
able  as  a  response  to  the  relative  magnitude  of  . 

Remark  4.3.  Theorem  4.1  is  directly  applicable  in  the  esti¬ 
mation  of  systems  reliability  when  two  or  more  components  are 


■1W  IJmJJ  l  W  UT JWlAflVW  LVWAW  **"  W  V-*  7-'T*  W 


•jOTiWW  rm  t\*  iv*  n» 


connected  in  series.  In  those  examples ,  often  our  problem  is  to 

estimate  I  6.  vhere  e.'s  are  scale  parameters  of  the  ith  life 
i 

distribution. 

Remark  4. 4.  Theorem  4.1  can  be  easily  extended  for  the  esti¬ 
mation  of  a  linear  transform  L68  where  L  ■  (  1  „  ),  an  r  *  p 
matrix. 


5.  NUMERICAL  STUDIES 

In  this  section,  we  study  PI  as  in  (3.8)  for  our  improved 
estimators  (2.2)  in  three  interesting  cases. 

In  the  first  two,  we  take  s.  *  1  and,  for  simplicity,  c.  *  c, 
£L  ■  -c/2,  i  ■  l,...,p.  This  implies  from  (2.3a)  that 

a  *  p  *(1  ♦  c/2)  and  the  improved  estimator  (2.2)  becomes 


6.(X)  -  a.X.  ♦  bx7c/2  H  x^1+c/2)/p 


,<1> 


1  i 

(2) 


(5.1) 


We  set  b  as  the  midpoint  of  this 


where  0  <  b  <  2p  d  /p-  d 
a  Za 

interval . 

Example  1.  F-d is t r ibut ion b .  Suppose  S  ■  (S  , ...,S  )  and 

2  1  P 

T  *  (T.,,,.,T  )  are  independent  where  S.  ^  o.  x_  and 

1  P  1  1  ^li 

2  2  •  2  2 
T.  ^  t.  y  »  i  *  l,...,p.  Our  problem  is  to  estimate  6.  *  o./i., 

x  x  ^  xxi 

i  *  l,...,p.  The  best  scale  invariant  estimator  of  6  ■  (6^,..., 6  ) 

is  6°(X)  •  (6?(X) . 6®(X))  where 

1  P 


6 . (X)  *  a .X.  ,  i  ■  1 , . . . ,p, 
l  ii 


(5.2) 


-1 


with  X.  *  S./T.  and  a.  *  («..  -  4)(n..  +2)  ,  i  ■  l,...,p,  and 

ill  l  n2i  li 


individually  the  X^  are  admissible. 

For  convenience  we  set  n^  ■  n^,  n^ 


n^  and 


i,o 


m  Q  n<  n«  n,  n^ 

•  (~>  Ti-±  *  a)H-~  -  o)/r(~)r(-~) 


(5.3) 


n. 
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A 


provided  -  n^/2  <  a  <  n^/2.  Thus  for  •  given  6,  we  cen  readily 
evaluate  (3.5),  (3.7),  and  ultimately  (3.8).  In  Table  I  we 
preaent  Pi  for  various  p,  c,  and  (n^,n^)  combinations.  6  was 
created  by  selecting  each  coordinate  randomly  within  the  given 
range.  The  results  for  a  typical  value  are  presented  along  with 
the  maximum  Pi  using  (3.9). 

Example  2.  Reciprocal  Beta.  Suppose  Y^  is  distributed  as 
n^*Be(2  +  e,l),  t  >  0,  i  *  1,2,. ...p,  and  we  seek  to  estimate 

n  *  m  (n,1,. . . ,n  *),  i.e.,  s  *  -1.  For  convenience  we  transform 
—1  P  —i 

to  X.  *  Y.  and  taking  0.  *  n.  we  have  at  0.  *  1,  f(X)  *  (2  +  e) 
ii  i  l  l 

•  X  The  best  invariant  estimator  of  0^  is  (1  +  e)  ^eX^ 

and  ia  admissible  from  the  discussion  at  the  beginning  of  Section  2 

since  e  >  0.  Using  (5.1)  and  noting  that  *  (e+2-ct)  (e+2),  a<  2, 

it  is  straightforward  to  calculate  (3.10)  which  becomes,  for 

c  >  -2, 


(i  ♦  c/2 r 
(1  ♦  c/2  ♦  eV 


(c  ♦  2)(c  ♦  2  ♦  c) 
(e  ♦  2  ♦  c/2)2 


(5.4) 


Using  (5.4)  we  see  that  for  p  large  if  either  e  or  c  “,  Tl  -*■  0 
while  if  c  and  c  are  close  to  0  nearly  1002  improvement  is  possible. 
Table  II  displays  results  for  small  to  moderate  p  using  (3.8)  with 
6  as  in  the  previous  example. 

Example  3.  Simultaneous  Exponential  and  Rectangular  Distribu¬ 
tions.  Here  we  consider  X^  arising  from  two  different  distribution¬ 
al  families.  Suppose  X^  have  p.d.f.  as  follows: 


f<xilei) 


,  -x./0. 

fi'1  »  1  1 

ei  e 


1 , .  •  • *Pi 


,  x.  >  0,  0.  >  0 

’  l  l 


p_+l,...,p  ,  0  <  x.  <  0.,  0.  >  0 

*1  *  ,r  *  l  ii 


(5.5) 


with  Pj  1,  P  i  2.  In  estimating  0  ■  (0j,...,0^),  the  best 

invariant  estimator  is  S^(X)  ■  (6?(X) . ^O(x))  where  6^(X)  “  a.X., 

1  P  ill 


«.»  •»*  I 


■<  (jk  h*  (.f  |, 


i.«.k .«i  »*4  »M  .4 


TABLE  I  (Continued) 


P 

>-* 

It 

o 

H 

N 

a 

o 

CM 

■ 

fH 

e 

nl  "  5 

nx  -  10 

p  -  2 

°2  '  5 

n2  “  10 

n2  -  20 

*»2  •  10 

n2  -  20 

(0,12) 

17.91 

23.14 

15.01 

24.55 

18.96 

6^  equal 

18.14 

23.44 

15.51 

24.87 

19.20 

P  *  5 

(0,12) 

15.78 

20.97 

15.82 

19.31 

19.01 

6^  equal 

34.87 

46.35 

34.97 

42.68 

42.01 

p  •  10 

(0,12) 

12.57 

19.84 

15.72 

15.60 

18.46 

6^  equal 

33.96 

53.59 

42.45 

42.13 

49.90 

P  -  30 

(0,12) 

10.07 

19.64 

16.12 

12.64 

18.59 

6.  equal 

29.53 

57.98 

47.60 

37.33 

54.90 

TABLE  II 
Reciprocal  Beta 

Percent  Relative  Improvement  in  Risk 


■  2 

,12) 

equal 

■  5 

,12) 

equal 


c  -  .01 

75.95 

76.20 


68.74 

91.74 


65.57 


56.36 

56.54 


56.00 

74.03 


54.43 


c  «  1 

10.21 

10.24 


14.43 

19.25 


15.28 


5* 

W 

I 

I 

1 

n 


I 

sv 


i' 
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In  order  to  insure  the  existence  of  necessary  expectations  in 
(2.2),  we  take  6  *  1/2  so  that  a  m  l/2p  whence  our  estimator  has 
ith  coordinate 


a.X.  ♦  bxV2  *xV2p 

li  l  j 


(5.6) 


In  (5.5)  “j  a  “  r(a  ♦  1).  i  Pj»  while  m^  Q  ■  (a  ♦  1)”  ,  i  >  p^. 
Again  taking  b  to  be  the  midpoint  of  (2.4),  we  can  calculate  6(6) 
exactly  using  (2.5).  Straightforwardly, 

P, 


R(e,{°)  ■  4  I  e”2  ♦  i  £  e.**2 
2  1  >  4  Pl*i  1 


so  that  Pi  in  (3.8)  may  be  obtained.  Some  values  are  presented 
in  Table  111  with  6  as  in  the  previous  examples.  Generally,  small 
improvement  is  observed.  In  fact,  if  p^  ^  2  and  p^  -  p^  ^  2,  one 
would  do  better  overall  by  creating  a  dominating  estimator  of  the 
form  (2.2)  for  the  first  p^  coordinates,  a  different  dominating 


estimator  of  the  form  (2.2)  for  the  remaining  p  -  p^  coordinates  and 


then  combining  the  two.  We  do  not  pursue  details  and  extensions 
here . 


TABLE  III 

Exponential /Uniform 
Percent  Relative  Improvement  in  Risk 


(0,12) 


pr  1 

CN 

N 

o. 

Pi"  2 

pr  8 

Pl*  5 

px«  io 

P  ■  2 

p  •  4 

tJ 

I 

I-* 

o 

p  -  10 

p  *  10 

p  -  20 

c  ■ 

0 

0.31 

1.46 

1.79 

2.53 

1.79 

2.71 

c  * 

1 

0.27 

1.26 

1.76 

2.39 

1.68 

2.65 

c  ■ 

2 

0.69 

1.09 

1.75 

2.32 

1.67 

2.55 

c  i 

-1 

0.69 

1.84 

2.22 

3.44 

2.96 

3.42 
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20.  ABSTRACT 


A  general  class  of  estimators  is  developed  for  improving  upon  best  scale 
invariant  estimators  of  two  or  more  arbitrary  scale  parameters  (or  powers  there¬ 
of)  for  arbitrary  positive  distributions  with  sufficient  moments  under  weighted 
squared  error  loss  function.  The  technique  is  to  compute  the  risk  difference  in 


terms  of  moments  of  the  distribution.  Some  conditions  are  obtained  under  which 


the  maximum  improvement  is  possible,  and  the  form  of  the  estimator  can  be  chosen 
to  achieve  this  maximum  along  any  specified  ray. 


The  result  is  then  extended  to  the  estimation  of  a  linear  transform  of  the 


parameter  vector.  Finally,  some  examples  are  given  with  numerical  calculations 
to  obtain  the  amount  of  risk  improvement . 
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